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Mapping the core mass function onto the stellar IMF: 
multiplicity matters 
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ABSTRACT 

Observations indicate that the central portions of the Present-Day Prestellar Core Mass Func- 
tion (hereafter CMF) and the Stellar Initial Mass Function (hereafter IMF) both have approx- 
imately log-normal shapes, but that the CMF is displaced to higher mass than the IMF by a 
factor F ~ 4 + 1 . This has lead to suggestions that the shape of the IMF is directly inherited 
from the shape of the CMF - and therefore, by implication, that there is a self-similar mapping 
from the CMF onto the IMF. If we assume a self-similar mapping, it follows (i) that F-N jr\, 
where r\ is the mean fraction of a core's mass that ends up in stars, and N is the mean number 
of stars spawned by a single core; and (ii) that the stars spawned by a single core must have an 
approximately log-normal distribution of relative masses, with universal standard deviation 
<x . Observations can be expected to deliver ever more accurate estimates of F, but this still 
leaves a degeneracy between rj and N Q ; and <x is also unconstrained by observation. Here we 
show that these parameters can be estimated by invoking binary statistics. Specifically, if (a) 
each core spawns one long-lived binary system, and (b) the probability that a star of mass M 
is part of this long-lived binary is proportional to M a , current observations of the binary fre- 
quency as a function of primary mass, £>(M,), and the distribution of mass ratios, p q , strongly 
favour 77 ~ 1.0 ± 0.3, N ~ 4.3 ± 0.4, cr ~ 0.3 ± 0.03 and a ~ 0.9 + 0.6; r\ > 1 just means 
that, between when its mass is measured and when it finishes spawning stars, a core accretes 
additional mass, for example from the filament in which it is embedded. If not all cores spawn 
a long-lived binary system, db/dM l <0, in strong disagreement with observation; conversely, 
if a core typically spawns more than one long-lived binary system, then N Q and rj have to be 
increased further. The mapping from CMF to IMF is not necessarily self-similar - there are 
many possible motivations for a non self-similar mapping - but if it is not, then the shape of 
the IMF cannot be inherited from the CMF. Given the limited observational constraints cur- 
rently available and the ability of a self-similar mapping to satisfy them, the possibility that 
the shape of the IMF is inherited from the CMF cannot be ruled out at this juncture. 

Key words: stars: formation - stars: mass function - stars: binaries - stars: statistics. 



1 INTRODUCTION 

Understanding the processes that determine the IMF, and why these 
processes appear to vary little with environment a nd metallicity, is 
one o f the main challenges in star formation (e.g. lElmegreen et alj 
2008). Recent observations of prestellar cores (i.e. the dense, grav- 
itationally bound condensations in molecular clouds that are pre- 
sumed to be destined to form individual stars or multiple systems) 
suggest that such cores have a mass function very similar in shape 
to the I MF, but shifted to higher masses by a factor of three to 



five (e.g.lMotte et al Jl 1998J: iTesti & Sargentll 19981: Ijohnstone et alj 



five (e.g. Motte et al. 1998; lesti & Sargent 1998; Johnstone et al. 
2000; Mott e et alj|200ll : Ijohnstone et alj|200ll : lstanke et alj|2006l : 
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Enoch et alj 12004 Ijohnstone & Baifvl 120061: I Young et alj |20Qi 

Nutter & Ward-Thompson 2007; Alves et al. 2007; Enoch et al. 



20081 : ISimpson et alj J2008I : iRathborne et al. 20ol; iKonvves et ail 

2010). The inference is that, in a statistical sense, there is a more-or- 



less self-similar mapping from prestellar cores to stars, and that the 
shape of the IMF is therefore simply inherited from the shape of the 
CMF. If true, this simply moves the problem to one of understand- 
ing the processes that determine the CMF, and why the outcome 
of these processes also varies little with environment and metallic- 
ity. In addition, we still need to understand how an individual core 
maps into an individual star or multiple system, and to what extent 
this process can really be viewed as statistically self-similar. 

The IM F has been evaluated bv lKroupal 1200 ll) and lChabriel 
J20031 2005). IChabrien finds that the IMF is well fitted with a log- 
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normal function merging into a power law at high masses. Theo- 
retical models and simulations of turbulent fragmentation suggest 
that the CMF may also approximate to a log -normal function merg- 
ing into a power law at high masses (Padoan & Nordlund 2002; 



IPadoan et alfeOOllHennebelle & Chabrieij|2008ll2009l) . 

However, these theories do not address the origins of stellar 
multiplicity. It is therefore timely to formulate the mapping be- 
tween core mass and star mass using simple distribution functions, 
so that the additional constraints imposed by stellar multiplicity can 
be taken into account. It turns out that these additional constraints 
can be accomodated quite easily, but strongly favour a mapping in 
which each core typically spawns N Q ~ 4 stars, with quite high 
efficiency, rj ~ 1 ; the individual stars spawned by a core have a log- 
normal mass distribution with standard deviation a Q ~0.3, and two 
of them end up in a long-lived binary system; the probability that a 
star with mass M ends up in a long-lived binary system is approxi- 
mately proportional to M. 

In the interests of simplicity, we ignore the high-mass power- 
law parts of the mass functions, and concentrate on the log-normal 
parts, since these are the parts that are best constrained by obser- 
vation, and they can be described with just two parameters: a log- 
arithmiqj mean and standard deviation. Therefore our conclusions 
are most pertinent to the mass range where this log-normal form ap- 
pears to be an acceptable approximation, say 0.03 to 3M . How- 
ever, it should be noted that our conclusions are not significantly 
changed if the high-mass power-law tail is included; this simply 
makes the maths more laborious and less precise. For a detailed 
discussion of the IMF and the eight parameters that m ay be needed 
to desc ribe it more completely, the reader is refered to Bastian et alj 
(2010). We limit our consideration of multiplicity statistics to (i) the 
binary frequency as a function of primary mass, and (ii) the distri- 
bution of mass ratios (for systems with Sun-like and M-dwarf pri- 
maries), again because these appear to be the multiplicity statistics 
that are most robustly constrained by observation. For the puipose 
of this paper brown dwarfs are counted as stars. 

In Section[2]we present the definitions and assumptions under- 
lying our model. In Section[3]we present the observational data we 
will use to estimate the model parameters. In Section|4]we describe 
the consequences of the model, using simple arguments; this dis- 
cussion pre-empts the results of the more rigorous statistical anal- 
ysis that follows. In Section|5]we describe how stellar statistics are 
evaluated for a particular model using Monte Carlo integration; and 
in Section [6] we define the parameter we use to measure the qual- 
ity of fit between a model and the observations. In Section [7] we 
describe the Markov Chain procedure for identifying the best-fit 
model parameters, and in Section [8] we present the results. In Sec- 
tion[9] we discuss the results and relate them to previous work, and 
in SectionfTOlwe summarise our main conclusions. 



2 THE MODEL 
2.1 Assumptions 

If one accepts tha t most stars are formed in cores (see e.g. 
iBressert et alj2O10h . the model has only four assumptions. 

Assumption I. The central portions of the CMF and the IMF 
are both log-normal. 



Throughout, all logarithms are to base 10. 



Assumption II. The mapping between them is statistically self- 
similar, which means that the distribution of the relative masses of 
the stars spawned by a single core must also be log-normal. 

Assumption III. When a core forms more than one star, two of 
these stars end up in a binary system that is sufficiently long-lived 
to contribute to the statistics of binaries in the field. All the rest 
ultimately end up as singles. 

Assumption IV. The relative probability that a star with mass 
M ends up in a long-lived binary system is proportional to M" . 

We note that these assumptions are not made because we be- 
lieve they are necessarily true, but because they are simple, and 
because it turns out that they suffice to fit all the observational con- 
straints that currently appear to be robust. 

In addition, we note that the long-lived binary systems that 
contribute to the field statistics are probably not the only ones that 
form in a core-cluster, simply the one s that survive its dissolution 
and subseq uent tidal perturbat i ons (e.g. Kroupal 19951) . There is evi- 
dence (e.g. lKohler et alj2008l ; IChen et alj|2013h that the multiplic- 
ity is much higher than in the field for young stars in some star 
formation regions, and includes a significant proportion of higher- 
order multiples. However, by the time stars arrive in the field, many 
of these systems are likely to have been destroyed, and the wider 
systems will continue to suffer attrition due to stochastic tidal per- 
turbations. 

In other words, there are two very different timescales in- 
volved in the mapping. The mean number of stars spawned by a 
core (A/q) and the mean total mass of the stars spawned by a core 
(hence the efficiency, rj ) are - ignoring stellar mass-loss, accre- 
tion, mass exchange and mergers - determined by processes that 
terminate once the core disperses, after at most a few mega-years. 
In contrast, the binary statistics are never completely settled. They 
evolve most rapidly during the birth throes of the core-cluster (the 
N stars formed from a single core) and during its dispersion, but 
they then evolve further due to interactions with other stars in the 
same large-scale cluster (here presumed to be an ensemble of stars 
formed from an ensemble of cores), and they continue to evolve, 
after the large-scale cluster dissolves, due to interactions with the 
ever changing background gravitational field (e.g. tidal perturba- 
tions from passing stars and molecular clouds). However these lat- 
ter perturbations are rare, and given that the typical field star has 
been in the field for many giga-years, its binary statistics should 
by now be well defined. Our model does not concern itself with 
the details of the dynamical evolution of the stars spawned by a 
single core; it simply focuses on the properties of systems that sur- 
vive to populate the field, posits that each core typically spawns 
just one such system, and shows that the observed binary statis- 
tics are reproduced well if this system tends to comprise two of the 
more massive stars spawned by the core. Other binary systems, and 
higher multiples, are spawned by a core, but we presume that they 
are disrupted on a timescale < 1 Gyr. One would expect the binary 
systems surviving in the field to be on average more massive and 
closer than the ones that have been disrupted. 



2.2 Input parameters 

Table [T] summarises the six model input parameters, viz. the log- 
arithmic mean, yu c , and standard deviation, <x c , of the CMF; the 
efficiency, rj, i.e. the fraction of a core's mass that is converted into 
stars; the mean number of stars, ,/V , spawned by a single core; the 
logarithmic standard deviation, cr , of the relative masses of the 
stars spawned by a core; and the dynamical biasing parameter, a. 
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Table 1. Input parameters regulating a single Monte Carlo integration, and the ranges of values admitted by the Markov chain. The prior for the Markov Chain 
is that, within these ranges, all values are equally probable. M c is the mass of a core, and \M„ | n ^j are the masses of the stars formed from a single core. 



Parameter 



He 



n 



Identity 



Minimum Maximum 



Arithmetic mean of log (M c /M ) 
Standard Deviation of log m \M c /M s J 
Mean Star Formation Efficiency m Core 
Mean Number of Stars Formed in Core 
Standard Deviation of log (M„/M ) 
Dynamical biasing parameter, d ln(p M )/d ln(M) 





-0.2 


+0.2 




0.3 


0.7 


*;f{M n }/M c 


0.0 


2.0 




1.0 


7.0 




0.0 


0.5 


W) 


-2.0 


5.0 



There are direct observational constraints on n c and cr c , but not, as 
yet, on rj, N , cr and a. 

We note that values of rj greater than unity are admissible, be- 
cause, between the time when the mass of a core is estimated and 
added to the CMF, and the time when its star formation is com- 
plete, the core can, and almost certainly does, grow in mass, for 
exa mple by accretion along the filament in which it is embedded 
(e.g. lSmith et al]|201lh . By the same token it is not necessary that 
all the stars spawned by a core form simultaneously. Indeed, numer- 
ical simulations suggest that some of the stars spawned by a core 
start to condense out of the filamentary material accreting onto the 
core, and may only reach the co re as it starts to disintegrate (e.g. 
lBatdl2012l ; lGirichidis et aT]|2012l) 

In addition, non-integer values of N are admissible. In such 
cases, we adopt the simple device of dividing cores between the 
integer values that bracket N Q . Thus, for example, N = 2.2 means 
that 80% of cores have N = 2 and 20% have N = 3. 

Apart from this device, we do not allow any variance in the in- 
put parameters, because to do so introduces extra input parameters, 
but does not significantly improve, or even alter, the fits obtained. 



2.3 Output parameters 

Given the four assumptions listed above, and values for the six in- 
put parameters, we can predict the IMF (which, being log-normal, 
is characterised by a logarithmic mean, fi s , and a logarithmic stan- 
dard deviation, <t s ), the binary frequency as a function of primary 
mass, b(M l ), and the distributions of mass ratio for systems with 
Sun-like and M-dwarf primaries, p 9 (M,). Our objective is to use 
observations of these output parameters (yU s , o~ s , b{M { ), p q (M { )) to 
constrain the model input parameters (fi c , cx c , r), N , <r ,a). 



Ijanson et al.1 J2012I): Bin 4, thos e of iRaghavan et alj l2010|): Bin 
5, thos e of lPreibisch et alj Jl999h ; and Bin 6, those of lMason et alj 
( 1998). For evaluating the quality of the fit, we give the first four 
bins equal weights, Wt. = 1/16, i = 1 to4, so that their combined 
weight is 1/4. The last two bins are given zero weight, because the 
stars in these bins are not strictly field stars] 2 ] Therefore these two 
bins should not influence the choice of best-fit model. They are in- 
cluded because - notwithstanding - the predictions of the best-fit 
model agree with them well (see Fig. [3}- 

For the distribution of mass ratios, q, we consider only 
primary-mass bins m = 3 and 4, sinc e these are the ones 
with relatively robust mass-ratio statistics ( Raghavan et alj 1201 01 : 
iReggiani & Meverl 1201 ll ; Ijanson et ai] l2012h . In both primary- 
mass bins, the d istribution of mass ratios appears to be flat 
( IReggiani & Meverll201ll) . We follow convention by allocating the 
mass-ratios to five equal bins, I = 1 to 5, so that bin I accom- 
modates va lues in the range .2(1- 1) < q < 0.2f. For primary- 
mass bin 3, Ijanson et alj ( 120121) conclude that, when allowance is 
made for selection effects, the distribution of mass ratios is flat, and 
therefore we simply set all the expectation values to V P3[ = 0.20, 



and all the uncertainties to U„ 



0.05. For primary-mass bin 



4, we adopt ex pectation values and Poisson uncertainties from 
IRaghavan et alj (2010). For all ten primary-mass/mass-ratio bins 
we allocate W Pm( = 1/40, so that their combined weight is 1/4. 



4 SIMPLE INFERENCES 

In Section[8]we present the results of a Markov Chain Monte Carlo 
analysis. Here we present simple arguments to preempt the main 
results of that analysis. 



3 OBSERVATIONAL DATA 

Table [2] summarises the expectation values, V x , uncertainties, Ux, 
and weights, Wx, accorded to the different observational param- 
eters, X, that the model seeks to predict. The weights determine 
the influence that different observed quantities exert on the overall 
quality of fit of a model (see Section[6]l, and by design they add up 
to unity. 

For the mean and stan dard deviation o f the IMF, // s and <r s , 
we use values informed by IChabrierl (2005), and since these two 
quantities appear to be quite well constrained by observation, we 
give them both a high weight, W^ = W„- =1/4. 

For the binary frequencies we consider six primary-ma ss bins. 
Binm= 1 (the lowest ma ss bin) represents t he res ults of lClose et al.1 
J2003h : Bin 2, those of iBasri & Reinersl J2006I) : Bin 3, those of 



4.1 The shift between the IMF and the CMF 

The mean mass of the stars that form from a given core are related 
to the mass of the core by the efficiency, r\ (the fraction of the core's 
mass that ends up in stars), divided by the number of stars formed 
from the core N . Hence the factor by which the peak of the CMF 
exceeds the peak of the IMF is given by 



lO^c^s' 






(1) 



If we adopt p. s = -0.6 ±0.05 (from |Chabrierll2003l). and fi n = 0.0 ± 
0.1 (from, e.g., lEnoch et~aT]|2006l ; [Young et al.ll2o6^ : lEnoch et ail 



2 The last two bins concern binaries with relatively high-m ass short-lived 
primaries in the Orion Nebula Cluster l lPreibisch et al.|l 999ft and a mixture 
of systems in clusters, associations and the field (Mason et al. 1998). 
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Table 2. Output parameters characterising the observed IMF and binary statistics (two lefthand columns), and parameters regulating the quality of the fit of a 
model to the observations (three righthand columns). Column 1 gives the name of the parameter in the model, and Column 2 its identity. Column 3 gives the 
observed Value (V) of this parameter, and Column 4 its Uncertainty (U). Column 5 gives the Weight (W) accorded to fitting the observed value. M s is the 
mass of a star from the whole ensemble of stars formed in a single Monte Carlo integration. The sources for the observational data are given in Section|5] 



Parameter Identity 



Observed Value 


Uncertainty 


Weight 


V »< = 


= -0.70 


u*- 


= 0.10 


v^- 


= 1/4 


*V 


= 0.55 


U <r S - 


= 0.05 


*v 


= 1/4 


Vb,= 


:0.20 


u b] = 


= 0.15 


w b] = 


: 1/16 


v h = 


:0.26 


u h = 


= 0.10 


w h = 


: 1/16 


v* = 


:0.34 


u h = 


= 0.04 


w h = 


: 1/16 


VH = 


:0.45 


u b4 = 


:0.03 


w h = 


: 1/16 


v b ,= 


:0.70 


u b ,= 


= 0.10 


w b5 = 


= 


v* = 


:0.85 


u bb ~- 


= 0.10 


w b6 = 


:0 


v», 


= 0.20 


"p3,l 


= 0.05 


Wp,, 


= 1/40 


'P3.2 


= 0.20 


Upv 


= 0.05 


w PX2 


= 1/40 


Vp„ 


= 0.20 


u P „ 


= 0.05 


W«.3 


= 1/40 


'P3A 


= 0.20 


"p3,4 


= 0.05 


w PXi 


= 1/40 


v P3 , 5 


= 0.20 


u Pv 


= 0.05 


Wp 3 , 5 


= 1/40 


v«,, 


= 0.10 


u m 


= 0.03 


w Pi , 


= 1/40 


v«, 


= 0.25 


u Pi , 


= 0.05 


W P 4, 


= 1/40 


v Pi , 


= 0.21 


u Pi ,, 


= 0.05 


»V 


= 1/40 


v Pii 


= 0.19 


u Pii 


= 0.04 


Wp4.4 


= 1/40 


Vp*< 


= 0.25 


U Pi s 


= 0.05 


w Pi , 


= 1/40 



h 

In 
In 
In 
In 
In 

P3,( 



in,i 



Mean of log (M s /M ) 

Standard Deviation of log (M s /M ) 

Multiplicity Frequency in (0.05, 0. 10) M Q 

Multiplicity Frequency in (0.05, 0. 17) M G 

Multiplicity Frequency in (0. 15, 0.60) M G 

Multiplicity Frequency in (0.8, 1.2) M G 

Multiplicity Frequency in (3, 50) M G 

Multiplicity Frequency in (20, 70) M Q 

Fraction of systems from primary-mass bin 3 in mass-ratio bin ( (( = 1 to 5) 



Fraction of systems from primary-mass bin 4 in mass-ratio bin C (C = 1 to 5) 



l2008l:lKonvves et al]|201Ch . we have F 22 4 ± 1 , whence 

N Q = Frj ^ (4±l)i/. (2) 

4.2 Raising the degeneracy between N and 77 

The degeneracy between N Q and r\ can be raised by considering the 
binary statistics. Two essential features of the binary statistics in the 
field are that - very roughly - the number of single-star systems is 
comparable with, but somewhat larger than, the number of binary 
systems, and the binary frequency is an increasing function of pri- 
mary mass (db/dM l > 0). The influence of these constraints can 
be understood with the following Gedankenexperiment. Suppose 
(purely for the sake of argument, and averaged over all masses) 
that 60% of systems are single and 40% are binary. This can be 
achieved in two ways. 

• A/q = 1 .4. In this case, 60% of cores have N = 1 and spawn 
singles, whilst 40% of cores have N = 2 and spawn binaries. This 
gives 0.26 < n < 0.44. However, it means that the components of 
binary systems are on average less massive than single stars, and 
therefore the binary fraction is a decreasing function of primary 
mass, which is the opposite of what is observed. 

• N = 3.5. In this case, each core spawns a binary system, 
but 50% have N = 3 so they spawn one extra single star, and the 
remaining 50% have N = 4 and therefore spawn two extra single 
stars. This gives 0.7 < n < 1. Moreover, provided a>0, the compo- 
nents of binary systems are now, on average, more massive than the 
single stars, and consequently the binary fraction is an increasing 
function of primary mass, as observed. 

There is therefore a strong preference for the larger value of 
AC , to ensure that db/dM^ >0. 



4.3 Standard deviation of the relative masses of the stars 
spawned by a single core 

Since the mapping of the CMF onto the IMF involves the convolu- 
tion of a log-normal CMF with a log-normal distribution of relative 
stellar masses, the logarithmic standard deviation of the IMF, cr s , is 
obtained by adding the logarithmic standard deviation of the CMF, 
<x c , and the logarithmic standard deviation of the relative stellar 
masses, cr , in quadrature, 

a- 2 = a 2 + a 2 . (3) 

SCO v ' 

A corollary of Eqn. ((3} is that - for a self-similar mapping - the 
logarithmic standard deviation of the IMF cannot be smaller than 
the logarithmic standard deviation of the CMF, 

cr s > cr c . (4) 

In interpreting this inequality, one must recognise that the log- 
normal CMF we are discussing here is one that represents a very 
large region embracing a representative ensemble of star formation 
regions; the log-normal CMFs inferred for individual star forma- 
tion regions can - and apparently do - have a range of means and 
logarithmic standard deviations, but together they cannot have a 
logarithmic standard deviation greater than that of the IMF and still 
admit a self-similar mapping. Since observations suggest <r c ~ cr s , 
this in turn implies that cr cannot be very large. 



4.4 Mass ratios 

Observations (Raghavan et al. 2010; Reggian i & Meverj uOllI ; 
Uanson et alj|2012 ) suggest that the distributions of mass ratio for 
binary systems having Sun-like and M-dwarf primaries are both 
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flat. In our model, this means firstly that cr can not be very small\j 
otherwise the range of stellar masses formed in a single core would 
be too narrow to produce low-q systems; and secondly that a can 
not be too large, otherwise the low-mass stars would have little 
chance of pairing up with the high-mass ones to prouce \ow-q sys- 
tems. 



If mass-ratio bin I of mass bin m contains C„,t systems, the corre- 
sponding mass-ratio probability is 

Pm,{ = -p~ ■ (10) 

The model can then be compared with the observational data. 



5 MONTE CARLO INTEGRATION 

For a single model (i.e. a fixed combination of the input parameters, 
/j c , cr c , 77, A/q, <t , a), we evaluate the stellar statistics as follows. 
First, a core mass, M c , is obtained by generating a Gaussian 
random deviate, Q, on (-00, +00), and setting 



M r 



10^c+£ CT c>M 



(5) 



Next, if A/q is non-integer, a value for N is obtained by generating 
a linear random deviate, £, on (0, 1), and putting 



N 



INT(A/- ), when £> N - INT(A/ ); 

INT(A^ G ) + 1 , when £ < N - INT(Af ). 



(6) 



Otherwise N = N Q . Then the masses of the N stars spawned by this 
core can be obtained by generating Gaussian random deviates, Q, 
on (-00, +00), and computing 



M. 



N 



lO^o 



(7) 



If N > 2, the integrated probability of each possible pairing of these 
stars (star n with star n') is computed, 



p , 

1 n.n 



1" £" {m°m°\ 

v=N-l v'=N . . 



(8) 



Finally, a linear random variate, £, on (0, 1), is generated, and the 
pairing whose integrated probability is just above £ is selected. 

This is repeated until a total of 10 7 stars has been created. 
Then the mean and standard deviation of the IMF, fi s and <t s , are 
computed (using the logarithms of the stellar masses). For each star 
that falls in one of the mass bins defined in Table [2] we note (i) 
whether it is the primary in a binary system, the secondary in a 
binary system, or a single star; and, if it is a primary, we also note 
which mass-ratio bin the binary falls in. If mass bin m contains P„, 
primaries and S m singles, the corresponding binary frequencjy is 



(9) 



*m "•" " n 



3 It turns out that finding a value of cr that is both small enough to satisfy 
Eqn. [5] and large enough to deliver low-g binaries, is the hardest constraint 

for the model to satisfy. 

We refer the reader to Reipurth & Zinnecker (1993) for a discussion of 
different measures of multiplicity and their various merits. The one defined 
in Eqn. J9} is in effect the multiplicity frequency, but we refer to it as the bi- 
nary frequency because we a re only considering binaries. As pointed out by 
Hubbe r & Whitworthl 420051) , the multiplicity frequency has the nice prop- 
erty that it is insensitive to whether a binary system is actually a higher- 
order multiple. We note parenthetically that there are in general other stars 
in each mass bin that are secondaries, but these do not explicitly affect the 
calculation of the b,„ . 



6 QUALITY OF FIT 

For each model (i.e. each Monte Carlo integration with a given set 
of input parameters, /j c , cr c , 77, N Q , cr , a), the quality of fit, Q, is 
given by a sum of terms, 

W x (X - Vx) 2 
AQ X = - "\ rl X> , (11) 



representing how well the model prediction for output parame- 
ter X (= fi s , <t s , b m [for/7! = 1, 2, 3, 4], p„,t [form = 3, 4; I = 
1, 2, 3, 4, 5]) matches with the observational constraints (see Ta- 
bleO. The overall quality of fit for a given model is then 



Q(fJ c ,ar c ,ri,N ,cr ,a) 
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u 2 

Pm,{ 



(12) 



The first two terms on the righthand side of Eqn. dl2t measure the 
ability of the model to reproduce the observed IMF (with an overall 
weighting of 50%); the third term (involving a single summation) 
measures its ability to reproduce the observed binary frequency as 
a function of primary mass (with an overall weighting of 25%); 
and the fourth term (involving a double summation) measures its 
ability to reproduce the distributions of mass ratio for systems hav- 
ing Sun-like and M-dwarf primaries (with an overall weighting of 
25%). A notionally perfect fit corresponds to Q = 0, and |<2| can 
be interpreted as the number of standard deviations by which the 
model departs from a perfect fit. 



7 MARKOV CHAIN 

7.1 Range of fi c and cr c 

HERSCHEL has allo wed much mo r e robu st evaluations of the 
CMF For example, iKonvves et all (2010) obtain (ji c ,cr c ) = 
(-0.22, 0.42) and (-0.05, 0.30) in - respectively - the entire Aquila 
field and the main Aquila subfield. Previously, E noch et alj (2006) 
have estima t ed (Hr ,CTr) = (-0.05 ±0.25, 0.50 ±0.10) in Perseus; 
lYoung et all d2006h have estimated (/^ ,oy) = 0.3 ±0.7, 0.5 ±0.4) 
in Ophiuchus; and lEnoch et alj {2008) have estimated (jU c ,cr c ) = 
0.00±0.04, 0.30±0.03) for an ensemble of cores from Perseus, Ser- 
pens and Ophiuchus. 

However, all these evaluations are convolved with a num- 
ber of uncertainties. In particular, the use of greybody fits to esti- 
mate mean dust temperatures, the mass opacity coefficients needed 
to convert fluxes into masses, and the distances assumed for the 
star-formation regions, all introduce uncertainty into the derived 
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Figure 1. The Q- values for all models along the Markov Chain that have Q > — 1, plotted against ji c ,cr c , r\, N Q , <r Q and a. 



masses, and hence into the yu c -values. cr c -values may be some- 
what less susceptible to these factors, but are affected by the fact 
that the cores on the low-mass side of the log-normal tend to be 
close to the completeness limit. Furthermore, we are here con- 
cerned with the values of yu c and cr c for the totality of all star 
forming cores, rather than those for a single region. 

To keep the problem tractable, we restrict the Markov Chain 
to values of fi c in the range -0.2 < /i c < +0.2. We discuss the 
consequences of taking fi c -values outside this range, in Section [8] 
For cr c we restrict the Markov Chain to values of a c in the range 
0.3 < cr c < 0.7. This choice is informed by the range of observation- 
ally inferred values, and by the fact that a c cannot exceed cr s . 



7.4 Markov Chain 

The ranges detailed above define the input parameter space, and our 
prior is that all values in these ranges are equally probable . The 
Markov Chain then starts at an arbitrary point in this space, and 
makes a biased random walk around the space. The components of 
a step are generated from Gaussian distributions. A step is always 
taken if A<3 = <2 NEW - <2 OLD > (i.e. if it results in an improvement 
to the fit). If A<3 < 0, the code generates a linear random deviate, 
_£, on (0, 1), and only takes the step if A(3 > ln(_£) (i.e. steps that 
produce a deterioration in the fit are less likely to be taken the larger 
the deterioration). The size of a step is scaled so that roughly half 
of all putative steps are not taken. 



7.2 Range of r\ and N Q 

We restrict the Markov Chain to values of n in the range < n < 2, 
and values of N in the range 1 < N Q < 7. Evidently, if a core ac- 
cretes very rapidly on the way to forming stars, higher n values are 
possible, but this turns out to be unlikely. The arguments presented 
in Section l4~2l suggest that higher values of N Q are inadmissible - 
unless each core spawns more than one long-lived binary, and the 
efficiency is increased still further (see Section|4T). 



7.3 Range of cr and a 

We restrict the Markov Chain to values of cr in the range < 
<x < 0.5, on the grounds that cr has to be smaller than <x s , and is 
probably also smaller than a c . 

We restrict the Markov Chain to values of a in the range -2< 
a < 5 . This choice is info rmed by numerical w ork on the dissolution 
of small-JV clusters (e.g. van Albada|[l9 68a b; McDon ald & Clarke! 
1 19931 : ISterzik & Durisenl Il998t iHubber & Whit worthl J2005L and 
references therein), which suggests that, if the dissolution of a core- 
cluster involves pure gravitational interaction between the stars, a 
single long-lived binary is the most likely outcome and it usually 
comprises the two most massive stars, which implies a s> 1. Con- 
versely, if there is dissipati on - for example, because the stars are 
attended by massive discs JMcDonald & Clarke! 19951) - other pair- 
ings become more likely, which implies a smaller a value. Flat 
mass-ratio distributions translate into a preference for small a. 



n c 


= -0.03 + 0.10 


°"c 


= 0.47 ± 0.04 , 


'/ 


= 1.01 ± 0.27, 


K, 


= 4.34 ± 0.43 , 


°"o 


= 0.30 ± 0.03 , 


<!- 


= 0.87 ± 0.64 , 


Q 


= -0.33, 



8 RESULTS 

From the Markov Chain, there is a single well defined Q peak in 
the parameter space explored, and the best fit is obtained with 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 

i.e. 0.33cr overall difference between the model and the observa- 
tions. 

The parameters of the CMF (yu c , cr c ) are compatible with those 
obtained from observation, although fj c has a rather large uncer- 
tainty, and we return to this point below. 

The effi ciency {rj) is mu ch higher than the values normally es- 
timated (e.g. lAlves et alj2007l) . 77 is also only j ust compatibleat the 
high end of the range calculated theoretically bv lMatzner & McKed 
(2000), but in their model these high values arise in cores that are 
intrinsically flattened (so that outflows can escape without sweep- 
ing up much core mass), rather than as a consequence of forming 
many stars. High notional efficiencies may be an indication that 
core s grow in mass wh ilst they collapse and fragment to form stars 
(e.g. lSmithetalj201ll) . 

The mean number of stars formed from a single core (N ) 
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Figure 2. Iso-<2 plots on principal planes through the best=rit model. On each row the ordinate (vertical axis) is the same model input parameter, from top to 
bottom in the order /i c , <x c , rj, N Q , cr and a. Along each row the abscissa (horizontal axis) cycles through the remaining model input parameters, in the same 
order. By scanning along a row one can see both which parameters are tightly constrained by the model, and which parameters are correlated. The false colour 
encodes the value of Q. (see bar on right of plot), and the contours correspond to £2+1, 2, 3, 4, 5. 
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is also higher than the values normally invoked. Mathematically 
this follows from the large 77 (see Eqn. Q3, but physically it also 
derives - inevitably, in a self-similar mapping - from the need to 
form binaries with a frequency that increases with primary mass 
(see discussion in Section l-Ot , 

The spread of stellar masses from a single core-cluster (o" = 
0.29 ± 0.07) is such that, if the stars are paired randomly, between 
33 and 56% of the resulting systems have mass ratio below 0.5. 
Thus, in order to produce a flat distribution of mass ratios, the dy- 
namical biasing parameter should not be too large, and this is what 
the model infers {a = 0.6 ± 1.0). 

In Fig.[T]we plot those values of Q generated along the Markov 
Chain that exceed - 1 (i.e. those models that deliver output parame- 
ters that are collectively within 1 a of the observations), against the 
different model input parameters. These plots show that the best-fit 
model input parameters are all well defined, apart from /j c . For- 
tunately fi c is already quite well constrained by observation, and 
likely to become better constrained in the future. If fi c were in- 
creased, the efficiency, 77, would have to be reduced proportionately 
(or each core would have to produce more than one long-lived bi- 
nary) — and vice versa. 

Fig. [2] illustrates how Q varies on planes through the best-fit 
solution, i.e. if just two of the model input parameters are varied. 
These plots are generated with a regular two-dimensional grid of 
models, and 10 7 stars per model. On each row, the ordinate is the 
same for all five plots, and the abscissa cycles through the remain- 
ing five input parameters. From the plots in the first row we see that 
yu c is weakly constrained, and also, from the second plot along this 
row, that if fi c is increased, 77 must be reduced proportionately. In 
all other cases, an horizontal scan of the plots in a row reveals that 
the parameter concerned (the ordinate) is very well, and uniquely, 
constrained by the observations. 

Fig. [3] presents the binary frequency as a function of primary 
mass, for the best-fit model, generated using 10 7 stars, along with 
the observational data used to constrain the model. We reiterate that 
we do not use the two higher-mass points, only the four lower-mass 
points. Notwithstanding, the model fits all six points well. 

Fig. [4] presents the mass ratio distributions for binaries hav- 
ing primaries in mass-bins 3 and 4. We see that there is acceptable 
agreement. The largest divergence occurs in the extreme bins. This 
is not surprising, given that, in the model, the components of a bi- 
nary system are drawn from a log-normal distribution of masses 
with a power-law weighting. 



9 DISCUSSION 
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Figure 3. The boxes represent the observational estimates of multiplicity 
frequency in different primary-mass intervals, as detailed in the text, and 
summarised in Table [2] The error bars represent the observational uncer- 
tainties. The dashed line shows the multiplicity frequency as a function of 
primary mass for the best-fit model. The unruly points at large M[ are due 
to small-number statistics. 
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Figure 4. The distribution of mass ratios for binaries having primaries in 
mass-bins 3 and 4. The plotted symbols with error bars represent the obser- 
vationally inferred expectation values and uncertainties: orange plus signs 
for mass-bin 3, and blue crosses for mass-bin 4. The histograms represent 
the model results: orange solid line for mass-bin 3, and dotted blue line for 
mass-bin 4. 



9.1 Critique of the model 

The critical assumption of the model is that each core spawns, on 
average, exactly one long-lived binary system, i.e. one binary sys- 
tem that survives to populate the field. If this assumption were re- 
laxed, in the sense that a core might spawn more than one long- 
lived binary system (say, on average S binary systems), then 77 and 
N would have to be increased (in proportion to S). Conversely, if 
not all cores were to spawn a binary system, 77 and N Q would have 
to be reduced, but it would then becomes impossible to reproduce 
the variation of binary frequency with primary mass, b(M [ ) — 
unless one were to introduce an additional parameter to allow the 
efficiency to be much higher for cores that spawn binaries than for 
those that don't. 

If the observed estimate for the overall binary frequency 



of low-mass field stars (i.e. binaries with primaries in the range 
(0.02, 2.0) M ) were to increase, this would reduce N , and conse- 
quently rj. For example, if the observed overall binary frequency of 
low-mass field stars were increased to 0.5, the model would require 
A/q-3 and 7/~0.8 ±0.2. 

It is difficult to see how the various standard deviation s could 
change much, unless <r s is very different from the IChabrierl ( 2005 ) 
value. If cr s were larger, <x c and cr could also be larger, and vice 
versa. 

If the distribution of mass ratios were skewed in favour of sys- 
tems with comparable mass, i.e. dp q /dq > 0, then cr would need 
to be reduced, and/or a increased (more dynamical biasing). 



© 2002 RAS, MNRAS OOO.fTliTol 



From the core mass function to the stellar IMF 9 



9.2 Previous theoretical work 

Some of the consequences of a self-similar mapping are explored 
by Clarke ( 1996), but with different distribution functions, and less 
emph asis on observational c onstraints . 

I Swift & William s (2008) develop a similar model to ours, but 
one which includes a power-la w extension to the CMF a t high 
masses, based on the analysis of iPadoan & Nordlunj 020020 . and 
which invokes somewhat different model parameters. They explore 
the consequences of varying the prescriptions for generating mul- 
tiple systems, and for sub-fragmentation of a core, but their work 
differs from ours in that they do not explore in depth the question 
of multiplicity and its variation with primary mass, and they do not 
draw any firm conclusions on the efficiency, or on the number of 
stars spawned by a single c ore . 

iGoodwin et alj (2008) explore the consequences of multiplic- 
ity for the mapping from the CMF into the IMF, and in particular 
the effect of multiplicity on the extremes of the IMF. Their pre- 
ferred model presumes that all cores spawn multiple systems, with 
the number of stars in a system increasing very slightly with the 
mass of the progenitor core (the model is therefore not strictly self- 
similar), and it has quite a low efficiency, r] = 0.27. They do not 
explore the issue of how such systems might subsequently evolve 
to produce singles, so they cannnot exploit the observed variation 
of bin ary fraction with primary mass. 

IGoodwin & Kouwenhovenl 12009) demonstrate that the map- 
ping from a log-normal CMF into an approximately log-normal 
System Mass Function (SMF) and from the SMF into an approxi- 
mately log-normal IMF admits a wide range of prescriptions for (i) 
how the efficiency, varies with the core mass (r](M c )), (ii) whether 
the probability that a core spawns a single or a binary depends on 
its mass (effectively N(M C )), (iii) and the distribution of mass ra- 
tios in such binaries. This concurs with our conclusion (see Section 
19.31 ) that, whilst there are many theoretical arguments for allowing 
the input parameters of the mapping to depend on the core mass 
(thereby rendering the mapping non self-similar), the effect on the 
IMF is so subtle that these dependencies cannot usefully be con- 
strained by the existing observations. 



of a simple free-fall argument, and Larson's scaling relations, this 
might be taken into the reckoning with^, = + 0.25. 

(ii) We have explored models in which the efficiency of star 
formation in a prestellar core depends on its mass according to 
rj oc Mp''. This is equivalent to including feedback from massive 
stars. Star formation is promoted by feedback from massive stars 
if Xn is positive, and suppressed if % n is negative. However, it is 
not known what the sense of feedback from massive stars is, on the 
scale of a single core. 

(iii) We have explored models in which the number of stars 
formed from a prestellar core depends on its mass according to 
N oc M C N . Negative values of^ w (a) skew the IMF towards high 
masses, and (b) increase the multiplicity frequency of high-mass 
stars and reduce the multiplicity frequency of low-mass stars. Pos- 
itive values of x N nave the opposite effects. It is difficult to believe 
that N does not increase with core mass (positive k,). However, 
this would completely undermine the original argument for a self- 
similar mapping between CMF and IMF, namely that the high-mass 
slopes of the CMF and IMF appear to be indistinguishable. More- 
over, in practice, the observational constraints can more easily ac- 
commodate the effects of negative x N - Either way, non-zero x N - 
values are not actually needed to fit the observational constraints 
we have invoked.. 

(iv) We have explored models in which the logarithmic range of 
stellar masses formed from a prestellar core depends on its mass 
according to <r oc M* <r . It is probably the case that only posi- 
tive values of Xir could be justified (i.e. higher-mass cores spawn- 
ing a greater logarithmic spread of stellar masses), but this is not 
needed to fit the observational constraints. Moreover, it suppresses 
the high-mass end of the IMF, which - in this purely log-normal 
model - is already too low. 

(v) We have explored the possibility that there is some variance 
in, for example, N Q , so that when N = 3 (say) not all cores spawn 
exactly three stars. However, firstly this introduces an additional 
model parameter, which should be avoided if possible, and sec- 
ondly it makes no significant difference to the results, unless the 
variance is extremely large, so we do not include it in the basic 
model. 



9.3 Additional model parameters 

We have considered the following refinements to the model. How- 
ever, none of them is justified, since none of them, either individ- 
ually or in combination, produces a significant improvement to the 
fit; in respect of items (ii), (iii) and (iv), Goodwin & Kouwenhoven 
( 2009) reached essentially the same conclusion, but on the basis of 
a very different model and less restrictive observational constraints. 
Necessarily, all these refinements would corrupt the self-similarity 
of the mapping. 

(i) We have explored models in which the lifetime of a prestel- 
lar core (i.e. the time during which a prestellar core is detected as 
such) depends on its mass according to t c oc M£'. Negative val- 
ues of Xt skew the IMF towards high masses, because low-mass 
cores are over-represented in the CMF. Conversely, positive val- 
ues of Xt skew the model IMF towards low masses, because high- 
mass co res are over-represented in the CMF. There is no consensus 
on this. lHatchell & Fulled |2008) have argued that more massive 
cores evolve faster than less massive ones, and are therefore under- 
represented in the CMF; this might be taken into the reckoning 
with^, = -0.25. Conversely. IClark, Klessen, & Bonnelll 120071) 
have argued that massive cores, being more diffuse have longer life- 
times, and are therefore over-represented in the CMF; on the basis 



We reiterate that we are not arguing that these additional ef- 
fects do not occur in nature. We are simply pointing out (a) that 
they are not justified by the currently available observational con- 
straints, that is, one can obtain a good fit to the observations without 
them; and (b) that they would corrupt a self-similar mapping. 



10 CONCLUSIONS 

We have developed a simple model to describe the mapping of the 
CMF onto the IMF. 

• The model has four assumptions: the central portions of the 
CMF and IMF are both log-normal; the mapping from the CMF 
onto the IMF is self-similar; if a core forms more than one star two 
of the stars end up in a long-lived binary; and the probability of a 
star of mass M being in this binary is proportional to M a . 

• The model has six input parameters: ji c and cr c are the loga- 
rithmic mean and standard deviation of the log-normal CMF; r] is 
the efficiency (i.e. the fraction of a core's mass that ends up in new 
stars); yV is the mean number of stars spawned by a single core; 
o" is the standard deviation of the log-normal distribution of rela- 
tive stellar masses spawned by a single core; and a is the dynamical 
biasing parameter. 
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• This model is able to fit the observed IMF, the observed bi- 
nary frequency as a function of primary mass, and the observed 
distributions of mass ratio for binaries having Sun-like and M- 
dwarf primaries. The best fit requires fi c = -0.03 ± 0.10, <x c = 
0.47±0.04, j]= 1.01+0.27, N =4.34±0.43, a Q =0.30±0.03, 
and a = 0.87 ± 0.64 . It fits the observations to within 0.25cr. 

We have not demonstrated, nor do we advocate, that the map- 
ping is necessarily self-similar, only that, if one assumes self- 
similarity, there is a simple mapping that fits the observational con- 
straints well and therefore - on the basis of Occam's Razor - should 
be given consideration. 

Moreover, if the mapping is not (at least, approximately) self- 
similar, then the notion that the shape of the IMF is inherited from 
the CMF must be abandoned. 

Either way, there is a question to be answered beyond under- 
staning the origin of the CMF: either why is the mapping self- 
similar, or why does the mapping, despite not being self-similar, 
produce an IMF with the same shape as the CMF? 

The self-similar model suggests that the efficiency of star for- 
mation within a prestellar core is significantly higher (r] n ^ 1.0 ± 
0.3) than has previously been proposed (e.g. t] ~ 0.3, lAlves et alj 
2007). It also suggests that most stars, including singles, are born 
in small groups of ~ 4. This contrasts with the conclusion of lLadal 
( 2006) that most stars , being single, are born in isolation. Interest- 
ingly |NakairruraetalJ ( 120121) have recently reported evidence that 
prestellar cores are more fragmentated than had previously been 
thought. If cores spaw n many stars, w e may see multiple outflows 
from some cores (e.g. IWu et al.| [2009), but these outflows do not 
have to disperse a large fraction of the core's initial mass, and can 
simply punch holes in the residual envelope. 
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